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ABSTRACT 

Helical structures are common in extragalactic jets. They are usually attributed in the literature to 
periodical phenomena in the source (e.g., precession). In this work, we use VLBI data of the radio-jet 
in the quasar S5 0836+710 and hypothesize that the ridge-line of helical jets like this corresponds 
to a pressure maximum in the jet and assume that the helically twisted pressure maximum is the 
result of a helical wave pattern. For our study, we use observations of the jet in S5 0836+710 at 
different frequencies and epochs. The results show that the structures observed are physical and not 
generated artificially by the observing arrays. Our hypothesis that the observed intensity ridge-line 
can correspond to a helically twisted pressure maximum is confirmed by our observational tests. This 
interpretation allows us to explain jet misalignment between parsec and kiloparsec scales when the 
viewing angle is small, and also brings us to the conclusion that high-frequency observations may show 
only a small region of the jet flow concentrated around the maximum pressure ridge-line observed at 
low frequencies. Our work provides a potential explanation for the apparent transversal superluminal 
speeds observed in several extragalactic jets by means of transversal shift of an apparent core position 
with time. 

Subject headings: galaxies: jets - hydrodynamics - instabilities - quasars: individual: S5 0836+710 



1. INTRODUCTION 

Jets in Active Galactic Nuclei (AGN) involve some of 
the most energetic processes in the Universe. The way 
they propagate and interact with the ambient medium 
and, in particular, their stability prop e rties have been 
thoroughly studied (see iHarded 120061 120111: iPeruchol 
[20T1I for reviews) . These jets are subject to the growth 
of different instabilities and many of the structures ob- 
served like knots, bendings, helices, have been inter- 
preted as a result of this physical process. Although 
there could be other ways to produce these structures, 
such as interactions with clumps of dense gas or pre- 
cession of the central engine, these in turn can also 
give rise to the growth of the instabilitie s via coupling 
to any of the unstable modes (see, e.g. iPerucho et al,| 
2p05;'Mizuno ct al. 2007; McKinncv & Blandford 2009; 
Mizuno et al.. ,2009: .Mignone et al., ,2010u ,Perucho et al.. 
20101 : IMizunoet al.ll201ir for numerical studies of cou- 
pling to Kelvin-Helmholtz and current-driven instabili- 
ties in relativistic flows). Our interest in understanding 
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how this process works in jets lies in the possibility of 
obtaining the physical parameters of the flow (velocity, 
gas density and s ound speed) , via modeling the observed 
struc tures (e.g., 'Lobanov fc Zensiii 120011 : iHardee et al.l 
I2005t 'Hardee & Eilek 2011). The Very Long Base- 
line Interferometry (VLBI) technique has provided high- 
resolution observations that produce highly detailed im- 
ages of the structur e of AGN jets, and that e ven resolves 
them transversally (|Lobanov fc Zensuil2001[ ). 

Radio images of many extragalactic jets from AGN 
present helical patterns on different spatial scales. There 
are a number of possible processes that can cause these 
structures, such as periodicity in the direction of ejec- 
tion of the flow due to precession. Precession can be 
produced by the gravitational effect of the accretion disk 
on the compact object in the center of the AGN or in 
a binary-black- hole object (see, e.g ., Lister ct al. 2003 
I Stirling et al. 2003; Lobanov fc Rolan d 2005; Bach et all 
l2006HSavoiainen et al, .2006) . If this periodical motion or 
any other initial perturbation of the relativistic jet cou- 
ples to an unstable mode of a magneto- h ydrodynamical 
instability like Kelvin-Helmholtz (KH, IPerucho et al.l 
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[200l lMizuno et alJ[2007l: IPeru cho et aII[20lol) or current 
driven (CD iMizuno eta l. 2009: 'McKinnev &: BlandfordI 
120091 iMignone et al.ll201Ql: iMizuno et al.ll201lD instabil- 
ities, the perturbation may grow in amplitude as it is 
advected downstream with the jet flow. Thus, the grow- 
ing helices can give us indirect information about their 
triggering process or eventually, about the physical pa- 
rameters of_the flow. 

iHarded ()2003[ ) compared the intrinsic and observed na- 
ture of helical structures via pseudo-synchrotron images 
of artifici ally bui l t heli cal jets from linear KH theory, 
following 'Harded (|2000( ). The resuhs of this work were 
applied to the jet in 3C 120 and allowed estimates to 
be made of some of the ph ysical parame t ers of the jet. 
This work was extended in iHardee et al.l ()2005[ ) , where 
the authors constrained the jet parameters more tightly 
and reproduced the basic observed structure of 3C 120 
by modeling the helical structures as KH modes. This 
typ e of analysis has also b een applied to the jet in M 87 
by IHardee &: Eilekl (|201lD who were able to constrain 
physical parameters in the jet and surrounding cocoon. 
Recently, iCaproni et al.l (|2009f ) have proposed a method 
(cross-entropy optimizer) to fit helical patterns produced 
by precession in jets. The parameters obtained by this 
method are the jet velocity, the jet viewing angle, and 
the jet opening angle. The authors generated artificial 
data-sets, which were successfully fitted by their method. 

Another problem related to observations of VLBI-jets 
is the measurement of supcrluminal transversal motion of 
components. iGongl (^2008.) has proposed non-ballistic mo- 
tion of physical components and misidentification as the 
origin of those measurements, although the latter is not 
a valid explanation for the case of densely time sampled 
jets ( e.g., the case of NRAO 150 and OJ287.'Agud o et al.l 
120071 2010, 2011). In this context, it is interesting to 
note that components fitted in radio-jets do not neces- 
sarily represent physical entities, but just a way to model 
the continuous structure of those jets. Only in some cases 
it is possible to associate the components that have been 
fitted to physical structures, such as perturbations cross- 
ing the jet or standing features, both usually related to 
shocks in the literature. 

In this paper, we present a deep analysis of multi- 
wavelength and multi-epoch observations of the helical 
jet in S5 0836-1-710 (0836+710 hereafter). Our analysis 
is based on the assumption that the maximum inten- 
sity of the radio emission coincides with the location of 
the pressure maximum in the jet. This is a reasonable 
assumption if the velocity of the emitting region of the 
fiow is fairly symmetric and there are no strong beam- 
ing asymmetries across the jet. The pseudo-synchrotron 
images produced in Hardee (2003), Hardee et al. (2005) 
and IHardee fc Eilekl (|2011( ) also indicate that the maxi- 
mum intensity in helically perturbed jets comes from the 
location where the pressure is at its maximum. Here, we 
present different tests that confirm this point from the 
observational data. Our analysis does not put any con- 
straint on the origin of the helices, but we hypothesize 
that the initial helical motion couples to either a KH or 
CD instability. 

The luminous quasar S5 0836-1-710 at a redshift z = 
2.16 hosts a powerful radio jet ex tending up to kilo- 
parsec scales ()Hummel et al.l I1992D . At this r edshift, 
Imas ~ 8.4 pc (see MOJAVE database and IWrighd 



2006 Q. VLBI monitoring of the source (jOtterbein et al.l 



1998f ) has yielded estimates of the bulk Lorentz factor 



7j = 12 and the viewing angle aj = 3° of the flow. The 
presence of an instability developing in the jet is sug- 
gested by the kink structures observed on milli arcsec- 
ond scales with grou nd VLBI (Krichbaum et al.l [1990.) . 
iLobanov et al.l ()1998f ) observed the source at 5 GHz with 
VSOFQ and also reported oscillations of the ridge-line. 
Identifying these structures with KH modes, they were 
able to derive an estimate of the basic parameters of the 
flow. High dynamic range VSOP and VLBA (Very Long 
Baseline Array of National Radio Astronomy Observa- 
tory, USA) observations of 0836-^710 at 1.6 GHz indi- 
cated the presence of an oscillation at a wavelength as 
long as ~ 100 mas (ILobanov et al.l [20061 ). which cannot 
be readily r econciled with the jet parameters giv en in 
ILobanov et all ([19981) . IPerucho fc Lobanovl (|2007D have 
shown that the presence of a shear layer allows fltting 
all the observed oscillation wavelengths within a single 
set of parameters, assuming that they are produced by 
KH instability growing in a cylindrical outflow. In this 
picture, the longest mode corresponds to a surface mode 
growing in the outer layers, whereas the shorter wave- 
lengths are identified with body modes developing in the 
inner radii of the jet. This result is reviewed here, in the 
light of additional information. 

The paper is structured as follows. In Section 2 we 
present the results from observations that will be used 
in later sections for our study. In Section 3 we describe 
the implications deduced from this set of multi- frequency 
and multi-epoch observations. In Section 4 we present a 
model to fit the ridge-line structures in relativistic jets. 
In Section 5 we discuss some properties of the jet in- 
ferred from the observations. Section 6 is devoted to the 
nature of the observed superluminal transversal motions 
in parsec-scale AGN jets in the context of our work. Fi- 
nally, we present our summary in Section 7. 

2. MULTI-EPOCH AND MULTI-FREQUENCY 
OBSERVATIONS OF THE JET IN 0836-1-710 

Observations of the jet at different frequencies and 
epochs are critical to our study. The different frequen- 
cies provide information on the stru ctures arising at dif- 
feren t spatial scales within the jet ([Perucho fc Lobanovl 
l2007f ). The different epochs provide a means to deter- 
mine structure motions. Table 1 shows the firequencies 
and epochs that have been us ed in this work: VL BA 
and VSOP at 1.6 and 5 GHz ( Lobanov et"al[ I2006D at 
three and two different epochs, respectively, two epochs 
at 1.6 GHz from EV1N0, one epoch of simultaneous global 
VLBI (including VLBA) observations at 2 and 8 GHz 
(01/1997, Pushkarev & Kovalev, submitted), two epochs 
from VLBA at 8 GHz, three epochs from VLBA at 22 and 
43 GHz, and 13 epochs, between 1998 and 2009, from the 

^ https;/ /www. physics. purdue.edu/astro/mojavc/, 

which uses ACDM Cosmology from WMAP 5 year 
results, and Ned Wrights Cosmology Calculator 
http: / /www. astro. ucla.edu/~wright/CosmoCalc. html 

^ VLBI Space Observatory Program, a Japanese-led 
space VLBI mission operated in 1997-2007 by the Insti- 
tute of Space and Astronautical Science, Sagamihara, Japan 
http: / /www. vsop.isas.jaxa.jp/top. html 

The European VLBI Network is a joint facility of Euro- 
pean, Chinese, South African and other radio astronomy institutes 
funded by their national research councils. 
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2 cm VLBA/MOJAVE database at 15 GHz. Figuresfflto 
[3] show several of these maps. 

The top two panels of Figure [T] show observations at 
1.6 GHz separated by about 10.17 yr but using two dif- 
ferent instruments. There is clear indication of a long 
wavelength oscillation and it can be seen that there are 
no significant changes at these large scales over this pe- 
riod. The bottom two panels show observations at 5 GHz 
separated by about 6.83 yr, both with the VLBA. There 
is some evidence for a long wavelength oscillation in both 
images but there are significant changes to the struc- 
ture with a suggestion of a shorter wavelength oscillation 
in the 2003 epoch at larger core distance (right bottom 
panel) . 

The top three panels of Figure [2] show observations at 
8 GHz, 19 and 16 months apart, spanning about 3 yr. 
The first image was obtained using Global- VLBI but the 
subsequent two images were obtained with the VLBA. 
These panels show only negligible difference, so there is 
little change in the jet structure over this 3 yr period. 
The bottom three panels show observations at 15 GHz, 
separated by 40 and 43 months and spanning about 7 yr. 
All three show evidence for some long wavelength oscilla- 
tion, longer than the observed jet, that does not change 
significantly over this period. The only significant differ- 
ence is that the first image clearly shows a shorter wave- 
length oscillation far from the core, whereas this shorter 
wavelength oscillation is weak at later times. 

The top two panels of Figure [3] show observations at 
22 GHz with 1.5 yr separation. Both show a straight 
jet along the first miliarcsecond, with some curvature 
appearing in the second image at larger distances from 
the core. The bottom two panels show observations at 
43 GHz with the same 1.5 yr separation and show a 
curved jet in the first epoch and a straighter structure in 
the second epoch. Note that the curved structure in the 
first epoch at 43 GHz could propagate out to larger scale 
with time and appear as the curved jet at larger scale in 
the 22 GHz second epoch. This would require a motion 
of about 0.2 mas in 1.5 yr (0.13 mas/yr ~ 3.5 c appar- 
ent velocity) . If the change in curvature were related to 
one quarter of a wave period, it would imply an observed 
periodicity of 6 yr. 

The ridge-line of the jet was obtained by determining 
the center of the jet emission (fitted by a Gaussian) at 
a given radial distance from the core. This calculation 
was performed at the different epochs and frequencies. It 
was done radially outwards to obtain a complete picture 
of the ridge-line of the jet. The computed ridge-lines are 
indicated as red diamonds on Figures [T][3] 

We have tested possible deviations in the computed 
ridge-line by obtaining it using two other approaches: 
from the location of the maximum emission and from the 
geometrical center of the emission profile in a transversal 
slice above a certain image rms cutoff level. As an ex- 
ample. Figure 0] shows the ridge- lines found by the three 
different methods for the 1.6 GHz image of the jet in 
1997 (also shown in Fig. [ij. The lines are overlaid for 
comparison and we find that there are no significant dif- 
ferences between the geometrical center of the emission 
profile, position of maximum emission and the center of 
jet emission as found by a Gaussian fit (as that shown in 
Fig. [l]). This coincidence will be reviewed below in the 
light of higher resolution observations at 15 GHz. 



In these comparisons, errors in the displacements are 
assumed to be one fifth of the FWHM of the resolving 
beam in the slice direction (this assumption is also used 
throughout the paper). We note that, in our specific 
case, this is a conservative assumption, as the positional 
errors are expected to be ~ 2FWHM/SNR (where SNR 
described the peak-to-rms signal to noise in the slice; c/., 
lCondonlfl997l) . SNR > 10 is obtained in all shoes cover- 
ing the entire relevant section of the jet, a conservative 
error estimate is justified to be used for the purpose of 
comparing the ridge lines obtained by the three differ- 
ent methods. With this approach, we find that all three 
methods are consistent with each other, hence implying 
that the resulting individual ridge lines are different by 
less than 20% of the FWHM. A more detailed measure 
of positional errors of the ridge-line, taking into account 
the effect of the SNR and the transverse width of the 
jet, will be employed in future work (Perucho et al., in 
preparation) in order to provide accurate comparisons 
with predictions from theoretical models. 

3. OBSERVATIONAL EVIDENCE 

In this section we discuss the ridge-line positions at 
different frequencies and epochs. In general, the appar- 
ent position of the opaque core depends on observing fre- 
quency and can shift with frequency on the sky (Lobano^ 
|1998[ ). However, we align images by apparent position of 
the core and do not take this effect into account while 
comparing the ridge-lines. This is justified by the small 
expected shifts (< 1 mas) as compared to the length of 
the structures that are compared (typically tens of mas). 
Se e for typical ain ount of the shift the statistical stud ies 
bv lKovalev eTall p008) and iSokolovskv et all (|20lll) . 

3.1. On the nature of ridge-lines 

Figures [S][7] display a combination of the radio emis- 
sion contours for different frequencies along with their 
ridge-lines for observations performed in 1997, 1998 and 
2003. Figure \5\ shows that the ridge-lines indicated at 
successively higher frequencies appear as oscillations on 
top of the ridge-lines obtained at the lower frequencies: 
the green 5 GHz and red 8 GHz ridge-lines wrap around 
the blue 2 GHz ridge- line. Differences between the 1.6 
and 2 GHz ridge-lines are within the errors. 

Figure [6] shows the same kind of behavior at higher 
frequencies in 1998, with the longest ridge- line belonging 
to the 8 GHz image and the ridge-lines at 15, 22 and 
43 GHz appearing as oscillations on top of the 8 GHz 
ridge-line. 

Finally, Figure [7] shows the broadest range of frequen- 
cies, from 1.6 to 43 GHz, with the bottom panels showing 
an enlarged image of the inner region. Again, the ridge- 
lines of the high-frequency images tend to lie on top of the 
lower frequency ridge-lines, except for the 15 GHz image 
ridge-line that deviates significantly from the 1.6 GHz 
ridge-line far from the core. However, at these distances, 
the signal-to-noise ratio is low in the 15 GHz image, 
which should translate into larger errors, and the differ- 
ence is always within the error of the 1.6 GHz ridge- line 
(~ 1 mas). 

In summary we find that the ridge-lines show remark- 
able coincidence independent of the epoch and the array. 
Significant differences would be expected if oscillation of 
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TABLE 1 

Observation Epochs and Frequencies 



YEAR 1.6 GHz 2 GHz 5 GHz 8 GHz 15 GHz 22 GHz 43 GHz 
1997 VLBA G-VLBI VLBA G-VLBI 



1998 VLBA ■■■ ••■ VLBA VLBA VLBA VLBA 

1999 ■ • ■ • ■ ■ • • ■ VLBA VLBA VLBA VLBA 

2000 ■ ■ ■ • ■ ■ • • ■ ■ • • VLBA 

2001 ■ • ■ • ■ ■ • • ■ ■ • • VLBA 

2002 ■ • ■ • ■ ■ • • ■ ■ • • VLBA 

2003 VLBA ■ ■ ■ VLBA ■ • ■ VLBA VLBA VLBA 

2004 ■ • ■ • ■ ■ • • ■ ■ • • VLBA 

2005 ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ • • VLBA 

2006 ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ VLBA 

2007 EVN • ■ ■ • • ■ ■ • • VLBA 

2008 EVN ■ ■ ■ ■ ■ ■ ■ ■ ■ VLBA 

2009 ■•■ •■■ •■■ ■■■ VLBA 



the ridge-line were array or epoch-dependent. This re- 
sult confirms the real nature of the structures observed 
in extragalactic jets by means of VLBI observations. 

If the jet is transversally resolved, the maximum in the 
emission, or the position of the ridge-line can be due to 
two main factors: 1) differential Doppler boosting caused 
by velocity gradients in the jet, or 2) an increase in the 
local pressure and number of emitting particles. If the 
position of this maximum is due to differential Doppler 
boosting, changes in the flow direction would have an in- 
fluence on the local emission from the different regions. 
This should result, depending on the internal velocity 
gradient, in irregular patterns, contrary to the results 
obtained. In addition, relatively large changes in speed 
and/or flow direction are required in order to get signif- 
icant differential Doppler boosting within the jet. Thus, 
the most reasonable hypothesis is that the ridge-line cor- 
responds to a pressure maximum. In the following sec- 
tions we present evidence that supports this interpreta- 
tion. 

3.2. Ridge-lines and transversal structure 

A straight axis can be defined for the jets, by a line 
starting at the core and defined by the symmetry axis of 
the the longest-observed wavelength oscillation (assumed 
to result from helical motion). Note that this line, which 
will be referred to as the jet axis or simply the axis, is not 
the same as the center or the maximum of the jet emis- 
sion at each radial distance, used above to determine the 
position of the ridge-line. Within our interpretation, the 
axis describes a mean or reference direction of the jet - 
the direction of flow propagation in absence of instabil- 
ities. This axis needs to be defined in order to be able 
to estimate the amplitudes and wavelength of individual 
instability modes, as done in Lobanov & Zcnsus (2001). 
The mean overall jet direction as defin ed by the jet 
axis, is modified by a helical instability (see lHardeell2000L 
for a discussion of the effect of different KH modes on 
jet structure and next sections). Variations of the ridge 
lines observed on smaller scales are likely to be caused 
by shorter unstable wavelengths. In order to estimate 
accurately wavelengths and amplitudes of these oscil- 
lations, the effect of the longest wavelength has to be 
taken into account, either by modeling it with an oscil- 
latory term or by assuming a mean local jet direction 
at each relevant spatial scale. As the small-scale ridge 
line variations occur on scales substantially shorter than 



the longest observed wavelength, we assume local mean 
jet directions for observations at different frequencies 
(relying on the fact that the spatial scales investigated 
depend largely on the observing frequency). If the jet 
axis is computed locally, for the different frequencies, the 
comparison shows that the absolute value of the jet posi- 
tion angle (x^,) increases with the frequency, where this 
angle is measured from the Declination axis, with 0° cor- 
responding to the positive values. 

The angle Xi^ does not change with time at each fre- 
quency. Figures [S] and IH] show this effect, with the po- 
sition angle of the jet ranging from 198° at 1.6 GHz to 
214° at 22 GH^, indicated by the dashed blue lines. The 
values of the different position angles at the different fre- 
quencies are listed in Table 2. 

Recall that there is a clear superposition of high fre- 
quency ridge-lines on top of the low frequency ridge-lines 
(Figs. [5][7]). This can be interpreted as the radio emit- 
ting regions at high-frequencies being concentrated in a 
relatively small region around the lower frequency ridge- 
line with the small scale oscillations seen at higher fre- 
quency developing on top of the long scale ones seen at 
lower frequency. This interpretation provides a natu- 
ral explanation of parsec-to kiloparsec scale jet misalign- 
ment for helically twisted jets pointing close to the line of 
sight. However, this raises a question as to whether the 
lowest-frequency emission observed is embedded in an 
even broader region, with lower or no emission. If this 
was the case, we would not be observing the whole cross- 
section of the jet at high-frequencies, but only a small 
region around the highest-emitting locations. Thus, the 
apparent helix would be tracing a high-emission region 
inside the jet, but not the whole jet. 

To test the possibility of a frequency dependent 
transversal structure in the emitting region, we com- 
puted the jet observed opening angles and widths for 
all the available frequencies at different pos itions along 
the jet axis following Pushk arev et al.l (j2009() . Note that, 
as stated above, the direction of the observed jet changes 
with the frequency and this should be taken into account 
to obtain these parameters at each frequency. Thus, the 
appropriate position angle of the frequency-dependent 
axis, Xi^j was used for each frequency. The results of this 
analysis are summarized in Table 2, where average values 

It is impossible to determine the jet position angle at 43 GHz 
due to its short length. 



Anatomy of helical relativistic jets: The case of S5 0836+710 



5 




Fig. 1.— Radio map of the jet in 0836+710 at 1.6 GHz (top) in 1997 and 2008 and 5 GHz (bottom) in 1997 and 2003. Isocontours range 
from 2 mjy to 1.45 Jy in the 1.6 GHz images and from 1.5 mjy to 1.08 Jy in the 5 GHz ones. The red diamonds indicate the position of 
the ridge-line. 



for each frequency are given. Not all epochs were used, 
as some were too noisy for this study. The errors come 
from the dispersion of the values for all epochs. We find 
that the opening angles are very similar at all frequen- 
cies, with an average value ofl3.8°±1.9. At the highest 
frequencies, 22 and 43 GHz, we have only considered one 
epoch (2003) due to problems with the other two epochs 
(1998, 1999). The resulting values are subject to larger 
errors as the opening angles are obtained with fewer cuts. 
If we disregard these high frequency results, the average 
opening angle is 12.1° ± 0.8, which at a 3° viewing angle 
implies an intrinsic opening angle of 0.63° ± 0.04. This 



value coincides with that obtained at 15 GHz, for which 
we have the largest sample and best quality data. The 
small or negligible differences between the opening angles 
at different frequencies tells us that either all transversal 
regions in the jet expand at the same rate or that, with 
the present resolution it is still difhcult to tell whether 
there is a transversal frequency profile for this jet. 

3.3. Possible secondary peak and implications 

Figure [TU] shows profiles of the jet brightness at dis- 
tances from the core between 6.4, 7.5 and 8.2 mas for 
different frequencies and epochs. These profiles do not 
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Fig. 2.— Same as[T]for 8 GHz (top) in 1997, 1998 and 1999 and 15 GHz (bottom) in 2002, 2005 and 2009. Isocontours range from 2.5 
mjy to 1.8 Jy in the 8 GHz images and from 0.5 mjy to 0.7 Jy in the 15 GHz ones. 



TABLE 2 
Observed Jet properties. 



1.6 GHz 5 GHz 8 GHz 15 GHz 22 GHz 43 GHz 



Average opening angle 

Position angle (x) 

Oscillation wavelengths (0-10 mas) 

Oscillation wavelengths (10-35 mas) 

Oscillation wavelengths (>35 mas) 



(12.3 ±1.2)° (13.5 ±1.3)° (12.1 ±0.8)° (10.5 ±0.9)° 18.1° 16.6° 

198° 202° 206° 210° 214° 

10 - 80 10 10 10 - - 

20 - 80 20 7-20 20 

40 - 80 - - - - - 



vary significantly with the position angle used, so we have 
taken xs, i-e., the position angle of the jet at 8 GHz as 
a characteristic value. The positions of these profiles are 
indicated in Fig. [T] In all the plots, the observed emis- 
sion maximum is displaced (by 1 — 2 mas) with respect 
to the reference axis, beyond the errors in the determi- 
nation of the peak. We will see in Section 1X5] that the 
measured values of the relative shift between epochs devi- 
ate from zero due to possible evolution of the ridge-lines. 
In most cases, single Gaussians provide good fits to the 
profiles and the center of jet emission coincides with the 
computed ridge-line. However, Figure [TU] shows that the 
transversal profiles at 8 and 15 GHz contain a secondary 



peak at the relative offset of —2 1 mas, which is not 

seen at the other frequencies at the same position. In 
order to verify that the second peak is real, a stacked 
imag(|3 from the MOJAVE database was used, and the 
same double peak structure in the profile was observed 
at this distance. In addition, if the image at 15 GHz 
is restored with the resolution of the 5 GHz data, the 
secondary peak disappears and the resulting profiles are 
very similar to the 5 GHz ones (Fig. [T0|) . Finally, the 
5 GHz image (shown in Fig. [Ij is accurately reproduced 
by the map of the jet at 15 GHz when convolved with 

^ This image is a composition from all the epochs, which rein- 
forces real features and smoothes noisy ones. 
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Fig. 3.— Same as [T] for 22 GHz (top) and 43 GHz (bottom) in 1998 and 1999. Isocontours range from 7 mjy to 5 Jy in the 22 GHz 
images and from 8 mJy to 5.8 Jy in the 43 GHz ones. 



the 5 GHz beam. 

At 15 GHz, where the two peaks are observed, the 
stronger peak is slightly displaced from the center of the 
emitting region (-3 to +4 mas, approximately), which 
is located around 2 mas to the right from the coordi- 
nate axis (jet axis, see Fig. [10]). The secondary peak is 
less important than the one to the right. This structure 
could be explained by superposition of a smaller ampli- 
tude elliptical mode and a dominant helical mode, or by 



limb-brightening of the jet. 

A growing instability in the linear regime implies a 
straight, expanding jet with the unstable pattern devel- 
oping in its interior or on the surface/la yer separating 
the jet from its surr ounding medium (^e.g.. lHardeell200Clt 
iPerucho et al.|[2"006[) . Thus, within this frame, the cen- 
ter of the jet emission and the ridge-line associated to 
the maximum emission do not necessarily coincide, if the 
latter is produced by the pressure maxima. The local di- 
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Fig. 4. — Radio map of the jet in 0836+710 at 1.6 GHz in 1997. The left panel shows the comparison between the ridge-lines computed 
from a Gaussian fit (red diamonds) and the geometrical center of the jet emission profile (blue diamonds). The right panel shows the 
comparison with the position of the maximum emission (black diamonds) . 




Relative R.A. (mas) Relative R.A. (mas) 



Fig. 5. — Radio maps at the 1.6 (black contours), 2 (blue contours), 5 (green contours) and 8 GHz (red contours) of the jet in 0836+710 
in 1997. The ridge-lines at the different frequencies are indicated by diamonds in the same colour as the contours and are shown without 
contours on the right. 



rection of the flow can change when the jet is affected by 
a hehcal mode and the maximum in pressure is displaced 
from the center of the flow, as show n in previous works 
dHardei Hop iPerucho et an[200l . The fact that the 
peak in emission at 15 GHz (Figure [TO)) is shifted from 
the center of the emitting jet supports the idea that the 
peak in emission does not coincide with the center of 



the jet flow and is more related to pressure differences 
in its cross-section. However, only if the jet is resolved, 
this feature is fully captured. This is the reason why, at 
1.6 GHz, where the jet is not well resolved, the different 
methods to determine the ridge-lines give similar results 
(see Section [2] and Figure |4]). 
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Fig. 6. — Radio maps at the 8 (black contours), 15 (blue contours), 
in 1998; the ridge-hnes at the different frequencies are indicated by 
contours on the right. 

3.4. Observed oscillation wavelengths 

Applying a rotation of the jet position angle to the co- 
ordinates of the ridge-line on the plane of the sky, the 
structure is brought to a horizontal axis (x-axis in carte- 
sian coordinates) along the positive Declination axis, 
making the detection of patterns much easier: 

X = Xob cos{2tt - x) - Hob sin{2TT - x) 

y^Xobsin{2TT-x) + VobCos[2TT-x), (1) 

where Xob and yob are the measured coordinates, x is the 
position angle of the jet in the plane of the sky defined in 
Section [3] and x and y are the rotated coordinates. Note 
that once rotated, the axis yob falls on the negative y axis. 
Figure [TT] shows the different coordinate systems used 
in this work. The coordinate system x, y is only used 
for wavelength detection purposes. In the next sections 
we will only refer to Xob, yob (observed positions) and 
Xh, yh, Zh (intrinsic coordinate system). 

Figure [iS] shows the transverse location of the ridge- 
lines relative to the jet axis (% — 198°) at a given epoch. 
The short ridge-lines fall on top of the long ones and 
do not show symmetry around the axis as they display, 
at most, half of the longest-observed wavelength, which 
determines the jet axis. In contrast. Figure [T51 shows the 
location of the ridge-line for different epochs at a given 
frequency (15 GHz), using the position angle of the jet 
at this frequency (xi5 = 210°) for the sake of clarity and 
to facilitate the detection of wavelengths in this plot. 

At 1.6 GHz, the ridge- line of the jet is dominated by a 
long wavelength oscillation of ~80 mas fFig. fT^ . A small 
scale modulation of 10 mas is also observed in the inner 
region (< 6 mas to the core, as indicated in the bottom 
right panel of Fig.[T2|). A ~ 20 mas structure is observed 
in the inner 40 mas at 2 and 5 GHz (e.g., left bottom 
panel in Fig.ll2p. Farther downstream, a high-amplitude 
(~4 mas) structure with wavelength ~40 mas is the most 
remarkable feature (e.g., upper right panel in Fig. I12p . 
There is an even shorter, small-amplitude wave, with 
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Relative R.A. (mas) 

22 (green contours) and 43 GHz (red contours) of the jet in 0836+710 
monds in the same colour as the contours and are shown without 

wavelength 5 mas, observed at 8 and 15 GHz in the inner 
15 mas, followed by another with wavelength 7-8 mas be- 
tween 20 and 30 mas (indicated in the bottom right panel 
of Fig. [T^ . Beyond 30 mas the data for the higher fre- 
quencies have low signal-to-noise. Figure [14] shows the 
ridge-lines at 1.6 GHz at different epochs (1997, 2003 
and 2008). Although the main structures are very simi- 
lar, the long wavelength is not so clearly observed in the 
last epochs and short-wavelengths seem to dominate the 
structure at larger distances from the core. 

The ridge-line at 1.6 GHz shows two peaks starting at 
80 mas from the core in the first epoch (Fig. [121 upper- 
left panel). The first peak does not appear in the next 
epochs (see Fig. [T4|) . In addition, there is a low signal- 
to-noise in this region. This leads to doubts as to the 
real nature of the first peak. The second increase in the 
amplitude (120-140 mas) is present in all three epochs. 

From the radio maps, we see that the jet emission 
drops when the ridge-line perturbation rises in the x-y 
plot at 1.6 GHz. This can be interpreted as the jet veer- 
ing slightly away from the line of sight. Moreover, the 
jet at higher frequencies is seen only up to 30-40 mas, 
i.e., exactly where this effect should be expected if the 
ridge-line veers away from the observer. This deviation 
is, nevertheless, very small because the amplitude of the 
perturbation is small compared to the axial wave-length 
{A ^ i mas and Aobs — 80 mas, see Fig. [T^. 

3.5. Wave speeds 

We have tried to measure the speed of the waves, both 
transversally and axially. To obtain a transversal motion 
we take fixed points in the Xob coordinate, defined by 
Equation [T] and Figure 111! and derive the displacements 
in the yob coordinate between two epochs. By doing this, 
we would expect an oscillation in the transversal veloc- 
ities if the apparent sinusoidal wave is moving spatially 
along the Xob-axis or if a standing sinusoidal wave pat- 
tern is growing temporally in amplitude. The first and 
second images in Figure [15] give the transverse velocities 
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Fig. 7. — Radio maps at the 1.6 (black contours), 5 (blue contours), 15 (green contours), 22 (red contours) and 43 GHz (violet contours) 
of the jet in 0836+710 in 2003. The ridge-lines at the different frequencies are indicated by diamonds in the same colour as the contours 
and are shown without contours on the right. The black lines in the top-left panel indicate the locations of the cuts shown in Fig. 1101 The 
bottom panels show an enlarged image of the inner region. 



that are found using the difference in ridge-Une position 
between 15 GHz VLBA observations made in 2002 and 
2003, and made in 2008 and 2009 (see Table 1). Fig- 
ure [15] shows transversal velocities from the radio core 
out to about 4 mas. The images for this part of the jet 
in 2002 and 2009 can be found in the bottom leftmost and 
rightmost panels in Figure[2j In FigurefTSlwe see that the 
transversal velocity is zero at the radio core since here it 
is assumed that the radio core is at a fixed position. In 
this regard note that the oscillation in transversal veloci- 
ties is not symmetric around zero and we will return to a 
more detailed discussion of this and other transversal ve- 
locity issues in Section |6l The oscillation pattern could 
be produced by a growing standing wave or a moving 
wave. However, the velocities obtained are larger than c 
and for a standing wave would imply amplitude growth 



at greater than light-speed and thus not possible. On 
the other hand, it is possible to find a transverse speed 
greater than c for a moving wave and we will show under 
what conditions in the next sections. Additionally, we 
observe an identical region of the jet that the oscillation 
pattern in transversal velocity has shifted between the 
epoch pairs along the by ^ 0.4 mas and includes 

only negative values for the transversal velocity. This 
shift indicates motion of the pattern along the Xo^-axis 
that would not occur for a temporally growing stand- 
ing wave. We will return to the reasons for the lack of 
similar maxima in the positive and negative velocities in 
Section [6l However, the errors in the ridge-line position 
are usually larger than the equivalent distance travelled 
by light in the interval between epochs. This results in 
values larger than c even if the true values are less than 
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Fig. 8.— 1.6 GHz map of the jet in 2003. The dashed-blue 
line indicates the position angle of the jet axis with respect to the 
positive Relative Declination. (198°). 
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Fig. 9.— 22 GHz map of the jet in 2003. The dashed-blue line in- 
dicates the position angle of the jet axis with respect to the positive 
Relative Declination. (214°). 

c. Actually, we have found that the maxima in transver- 



sal velocity are larger for shorter time-intervals between 
epochs, where less real motion is expected and differences 
in positions are dominated by errors. Thus, we cannot 
conclusively rule out a temporally growing standing wave 
based on the transverse velocities in Figure 1151 but ex- 
amination of the differences between the 2002 and 2009 
15 GHz panels in Figure [5] allows us to conclude that 
a moving wave pattern is the most likely cause of the 
transversal velocities. 

It is also possible to compute the axial speed, using 
fixed Hob position and computing the difference in the 
projected axial coordinate between epochs. We would 
expect a constant value for this speed if the wave is mov- 
ing. In fact we do not obtain a constant velocity. Al- 
though the evolution of the ridge-lines at 15 GHz (see 
Fig. [T3| shows that the small peak at 4 mas in the left 
top panel seems to move to the right with time, no regu- 
lar velocity pattern has been found in relation with this 
motion. The reason for this might be errors in determin- 
ing the axial locations due to blending of the axial and 
the second transversal spatial coordinates on the plane 
of the sky. In addition, whereas in the case of only one 
dominant wavelength, the measure could be easy, the 
superposition of the various different wavelengths with 
different potential speeds would lead to variability in dis- 
placements from epoch to epoch. In the next sections, we 
provide the equations that couple the transverse and ax- 
ial speeds with the amplitude and wavelengths. However, 
a detailed derivation of these equations will be presented 
elsewhere (Perucho et al., in preparation). 

4. FROM OBSERVATIONS TO THEORY 
4.1. Observed and intrinsic structures 

Past work has shown that oscillations like that 
observed in 0836+710 can be produced by three- 
dimensional helices developing in a relativistic outflow. 
In this case the position of the ridge-line will trace the 
accompanying helically twisted pressure maximum along 
the jet, e.g., the evidence presented in Section |3l In 
order to model these structures we have to take into ac- 
count projection and relativistic effects due to the possi- 
ble propagation of the waves with a high velocity. 

From Figure [11] we can see that the coordinate yob 
is not affected by projection or light travel time-delay 
and coincides with the corresponding coordinate of the 
intrinsic helix yh- However, coordinate Xob is affected by 
both projection and light travel-time effects and depends 
on both the transversal coordinate Xh and on the position 
Zfi along the axis of the helix and the jet itself. 

The observed helical wavelength is related to the in- 
trinsic on43 via the viewing angle, a, and the intrinsic 
wave speed, following 



Aof) 



sm a 



(2) 



(1 — cos a) ' 
The intrinsic equations describing a helix moving along 



^ We will call intrinsic a perpendicular position (viewing angle 
of 90°) of an observer at rest with respect to the jet flow. 

^ cf. Perucho &: Lobanovl |2007), where cosmological effects were 
mistakenly considered: the redshift effect in the frequency also 
affects the measured velocities, but not the measured wavelengths. 
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Fig. 10.— Profiles and Gaussian fits for 1.6 and 5 GHz (2003 epochs), 8 GHz (1998) and 15 GHz (2006), from left to right. The location 
of the cuts is shown in the top-left panel of Fig. [7| (the first two panels here are taken from the epoch in that Figure, whereas the other 
two have been plotted together for the sake of clarity). In the left panel (1.6 GHz), the cut is done 7.50 mas from the core along the axis 
(PA=— 153.4°), the peak intensity is 482.5 mjy/bm, the peak relative position is 0.591 mas, and the FWHM of the Gaussian is 6.630 mas. 
In the second panel from the left (5 GHz), the cut is done 7.40 mas from the core along the axis (PA=— 153.4°), the peak intensity is 
45.5 mJy/bm, the peak relative position is 1.454 mas, and the FWHM of the Gaussian is 2.533 mas. In the third panel (8 GHz), the cut is 
done 6.40 mas from the core along the axis (PA=— 153.4°), the peak intensity is 7.3 mJy/bm, the peak relative position is 1.298 mas, and 
the FWHM of the Gaussian is 2.362 mas. In the right panel (15 GHz), the cut is done 8.20 mas from the core along the axis (PA=— 153.4°), 
the peak intensity is 2.6 mJy/bm, the primary peak relative position is 1.859 mas, and the FWHM of the Gaussian is 1.31 mas. 

jected onto the sky plane {xob,yob) where the helical axis 
coordinate, Zh, and Xh are to be projected onto the sky 
plane x-axis and the helical axis makes an angle, a, with 
respect to the line of sight (which lies in the Xh, Zh plane). 
The relation between the observed {xob, Hob) and intrinsic 
{xhiVh, Zh) coordinates is the following (Perucho et al., 
in preparation): 




Xob = Xh{zh) 



1 — COS a 



Zh- 



sma 



' {1 — Pw COS a) 

Vob = Vh ■ (4) 
which, in the absence of motion reduces to: 



Xob — Xh cos a + Zh sm a 
Vob = Vh, 



(5) 



observer 



Fig. 11. — Coordinate systems used in the text. The intrinsic 
coordinates of the helix are x^, Vh and z^- The observed coordi- 
nates are x^t and y^i,. x, y, x^b, and the axis including y^t and y/^ 
fall on the plane of the sky. Note that this image is mirrored with 
respect to the radio maps in the paper. 



the Zh axis with speed Vw are: 

Xh = Ao e^"^^' cos ( + ipo 

Vh = Ao e^"/-^' sm ( '- + 



(3) 



where is the initial amplitude of the helix, Ai is the 
growth length of the helix if it couples to an instability, 
A is the wavelength, and <po is the phase. These four 
variables plus the wave velocity are the unknowns that 
determine the structure of the wave. Finally, t is the 
time, and Zh is the axial coordinate. This helix is pro- 



The observed ridge-line points could then be fitted to 
a helix by means of a minimizing method, using the 
above equations. However, Equation (j4]) embeds a key 
degeneracy in the problem: the observed wavelength is a 
combination of the wave velocity and the intrinsic wave- 
length. Thus, there are different combinations of both 
parameters that may result, for a given viewing angle, 
in the observed wavelength. In order to uniquely iden- 
tify the intrinsic wavelengths it is necessary to measure 
the wave speed of the different structures. A simple way 
to do this would be to measure motions in the ridge- 
line between epochs by following visible features such as 
maxima or minima and relate those to intrinsic motions. 
However, we have seen in Section 13.51 how difficult this 
may be with the present accuracy of the VLBI observa- 
tions. Thus, we can only find limits to the parameters 
implied by the structure. 

An additional complication arises when more than one 
periodical structure is found in a jet. In this case each 
periodical structure has to be independently fitted if 
it propagates with different speed. The reason is that 
the relativistic and projection effects, which depend on 
the propagation velocity, determine the position of the 
displacements along the observed axis following Equa- 
tion (j4]) , which makes independent parametrization nec- 



Anatomy of helical relativistic jets: The case of S5 0836+710 



13 



E 



crosses: 1.6 GHz 
triangles: 2.3 GHz 




+ 

+ 
+ 




1 


+ 


+ 
+ 






+ 

+ J 




t** 80 mas 










+, 





50 100 
X (mas) 



150 



U5 



crosses: 1.6 GHz 
diamonds: 5 GHz 



2 - 20 mas 




40 mas 



50 100 
X (mas) 



150 



triangles: 2.3 GHz 
diamonds: 5 GHz 




triangles: 5 GHz 
diamonds: 8 GHz 




Fig. 12. — Transverse displacement of the jet ridge-line around the jet axis at several frequencies in 1997. 1.6, 5 and 8 GHz images were 
obtained with VLBA data, whereas the 2 GHz one was obtained with Global VLBI. We indicate the wavelengths by the observed sizes 
between maxima or minima (half a wavelength in the case of the 10 and 18 mas oscillations). The black or red crosses show the error 
bars (one fifth of the beam size in the direction of the displacements), at randomly selected points, to be applied to each data point in the 
ridge-line at a given frequency. The errors are not shown at all points for the sake of clarity. 



essary. Once the velocity and intrinsic wavelength of each 
wave are known, their contributions to the ridge-line dis- 
placement, Uob in Equation m at a given location, Xob, 
have to be added and then compared to the observed 
ones. 

A more detailed description of the modeling of rela- 
tivistic helices will be published elsewhere (Perucho et 
al., in preparation). 

5. JET PROPERTIES DERIVED FROM THE 
OBSERVATIONS 

5.1. The 40 mas wave 

In all three ridge- lines observed at 1.6 GHz, a 40 mas 
wavelength between a; ~ 80 — 120 mas (see Fig. [U, in- 
dicates a shortening of the wavelength dominating the 
shape of the ridge-line and deserves some discussion. 
This structure can be explained in one of the following 
three ways: a) the amplitude of a shorter wavelength in- 
creases to larger than that of a longer one, b) the angle 
to the line of sight decreases and the long wavelength 
is observed with a shorter wavelength, or c) the longer 
wavelength is decelerated. Options b and c follow from 
Equation [2] that relates the observed and intrinsic wave- 
lengths, taking into account projection and relativistic 



effects. 

We discuss here the different possibilities 

• a: The longer initial wavelength could be due to 
a forced change in the jet direction that due to its 
low frequency couples to a long wavelength slowly 
growing KH mode. The shorter wavelength is the 
result of a faster growing, possibly resonant, KH 
mode, whose amplitude grows sufficiently rapidly 
to become larger than that of the forced oscillation. 
Thus, this is a viable option. 



b: We can find the required decrease in the an- 
gle to the line of sight for the 80 mas wavelength 
to change to the 40 mas wavelength using Equa- 
tion [21 and keeping the intrinsic wavelength and 
wave speed constant. For small viewing angles 
Equation[2]becomes 1— /J^uj^i cos ai ~ 1— /3u),2 cos Q!2 



* Considering that the changes discussed occur already at large 
de-projected distances from the core and that there are evidences 
of KH ins tability to dominate at distances lar ger than a few parsecs 
(see, e.g.. IHardei [200^ . [20Tll : IPeruchcll20111 '). we will restrict the 
discussion to this kind of instability. 
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Fig. 13. — Transverse displacement of the jet ridge-line at 15 GHz and different epochs around the position angle of the jet at this 
frequency, ranging from 1998 to 2009. In the upper left panel, black crosses represent the ridge-line from March 1998, red triangles from 
June 2000 and blue diamonds from June 2002. In the upper right panel, black crosses represent the ridge-line from March 2003, red triangles 
from September 2004 and blue diamonds from May 2005. In the lower left panel, black crosses represent the ridge-line from October 2005, 
red triangles from September 2006 and blue diamonds from Augus t 20 07. In the lower right panel, black crosses represent the ridge-line 
from July 2008 epoch and red triangles from May 2009. As in Fig. 1121 the blue crosses show the error bars (one fifth of the beam size in 
the direction of the displacements), at randomly selected points, to be applied to each data point in the ridge-lines. The errors are not 
shown at all points for the sake of clarity. 



with an error of < 10%, and: 



sm ai 
sina2 



(6) 



This means that a2 — ai/2, and a decrease in the 
angle to the hne of sight down to one half of its orig- 
inal value is needed. This could be produced by a 
change in the conditions of the ambient medium 
when the jet exits the ISM of the galaxy. However, 
this would imply the presence of a standing shock 
where the jet changes direction. Such a shock is 
not observed and there is no hint that the jet di- 
rection is appreciably changed. Thus, this option 
is unlikely. 

c: A decrease in the velocity of propagation of the 
80 mas wave could also lead to the observed wave- 
length decrease. It is well-known that the pertur- 
bations propagating in a flow cannot have velocities 
larger than that of the underlying fluid. For exam- 
ple, if the flow decelerates, the wave would slow as 
well. Now using Equation [2] and keeping the intrin- 



sic wavelength and the viewing angle constant we 
have: 



(1 ~ (3w.2 COS a) _ Xobs.i 
(1 - /3tu,i cosa) Xobsa 



(7) 



which should be solved for l3w.2- For example, 
taking i — 0.9 for the inner region, a — 3°, 
Ao&s,i = 80 mas and Xobs.2 = 40 mas, we ob- 
tain f3w,2 — 0.8, whereas for /S^^i — 0.6 we get 
/3tu,2 = 0.2. This option allows us to place con- 
straints on usingAo6s,i(l-/3„,iCosa) < Xobs,2, 
to obtain (3^,1 > 0.5. It is known that a flow can be 
decelerated when a helical perturbation grows to 
nonlinear values and disrupts the collimated flow 
resulting in c onversion of kinetic i nto th ermal en- 
ergy (see e.g.. lPerucho et al.ll2005ll2010D . The un- 
coUimated structure of the 0836-1-710 jet at kilo- 
parsec scales indicates that this is a viable option. 

Figure [H] showed that the long wavelength is not 
clearly observed in the later epochs at 1.6 GHz and short- 
wavelengths seem to dominate the structure at larger dis- 
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Fig. 14. — Rotated ridge-lines obtained for 1.6 GHz in 1997 (top), 
2003 (center) and 2008 (bottom). As in Fig. [12] the red crosses 
show the error bars (one fifth of the beam size in the direction of 
the displacements), at randomly selected points, to be applied to 
each data point in the ridge-lines. The errors are not represented 
for all points for the sake of clarity. 



tances from the core at the later epochs. It seems likely 
that some combination of options a and c are favored by 
the observations. 



5.2. On the presence of a shear-layer in the jet 

IPerucho fc Lobanovl (j2007f) presented results that 
pointed towards the jet in 0836+710 having a signifi- 
cant velocity shear layer. For this study, the authors 
used as a set of j et parameters those r esulting from a 
previous work by iLobanov et al.l ()1998[ ): Lorentz fac- 
tor 7 = 11, Mach number Mj = 6, and jet/ambient 
density ratio Pj/pa = 0.04. They studied the different 
characteristic wavelengths growing in a jet with a thin 
shear-layer (~ 10% of the jet radius), and a thick one (~ 
60% of the jet radius). The result showed that the ob- 
served wavelengths in the jet structure could be better 
explained transversally stratifie d jet. A m istake in the 
equation used in that work (Pe rucho fc L obanov 201 1]) 
led to misidentification of the modes in the original work, 
but the conclusions were shown to remain valid. 

In the present work, we have been able to accurately 
measure the wavelength of the longest mode to be closer 
to 80 mas than to 100 mas. A correct treatment re- 
duces the values of the wav elengths obtained in Table 1 of 
IPerucho fc Lobanovl (|2007f ) by a factor of ~ 3. This im- 
plies that their solutions both with and without a shear 
layer cannot explain the observed structures, r eported to 
be 100 mas and 7.7 mas (|Lobanov et al.lll998D . with the 
jet parameters used to obtain the results shown in Fig- 
ure 1 of Perucho fc Lobanov (2007). We reproduce here 
the solutions for the sheared jet in Fig. [T6l 

Taking 80 mas for the wavelength, 3° for the viewing 
angle, and that the velocity of this long wave has to be 
/3w > 0.5 (see Section l5.ip . we obtain an upper limit to 
the intrinsic wavelength using Equation [2] of 



ini,80 mas 



< 



1 - 0.5 COS 3= 
sin 3° 



80 mas = 765 mas. (8) 



In terms of the solutions presented in 
IPerucho fc Lobanovl (|2007|). this is A,>,f sn „ip.; < 45i?i, 
II ■ 



with Rj = 17 mas ([Perucho fc Lobanovl l2007'). Now the 
maximum growth rate of the surface mode occurs at 
160 Rj for the sheared jet and at 105 Rj for the vortex 
sheet jet, which is clearly different from the obtained 
intrinsic wavelength. For the sheared jet, the maximum 
growth rate of the first body mode is at the frequency 
uj ~ 0.1c/ Rj and A ~ 42i?j. Taking into account that 
its velocity is ~ ~ 0.5 c (computed as ~ Re{u]/k)), 
we obtain a value very close to Amt somas- We note that 
the maximum growth rate of the first body mode for the 
vortex sheet solution is at A ~ 21i?j, with ~ 0.6 c, 
farther away from the required properties to explain 
the observed wavelength. While the first body mode 
is not the fastest growing one in the solution, it could 
be triggered and generate large deformations in the jet, 
which a re not expected to be caused by high-order body 
modes (jPerucho et al.ll2005l l2007l l2oTol ). We note that 
this interpretation is very dependent on the measured 
jet radius. 

At the same triggering frequency (w ~ 0.\c/Rj) in 
the sheared case, high-order body modes could produce 
the 7-8 mas wavelen gths (identified with the 7 .7 mas 
wavelengths studied in IPerucho fc Lobanovll200"7[ ). This 
wavelength could really be ~10 mas, resulting an intrin- 
sic A^yii^iOmas 

~ 9i?j, obtained using 3° viewing angle 
and ~ 0.2 c (~ Re{uj/k)). The fact that both growing 
modes would be triggered at the same frequency sim- 
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Fig. 15. — Transversal velocities measured from 15 GHz VLBA data between epochs 2002-2003 and 2008-2009, respectively. 




Fig. 16. — Solutions of the linear problem of KH instability for a sheared jet (see lPerucho fc Lobanovll20071 . for details). The left panel 
shows the real (upper curves) and imaginary (lower curves) parts wavenumber with wave-frequency, whereas the right panel shows the 
wave velocity (R,e{Lj / k)) with wave- frequency. 



plifies the nature of the driving mechanism or mecha- 
nisms to a single one. In this case, the small wavelengths 
can be triggered as h armonics of the longest excited one 
(|Perucho et all 120051 ) . The new resulting frequency is 
a factor of 4 larger than that given in PL07 for the 
longest wavelength. Thus, the driving period would be 
now Trfr — 1-4 X lO"^ yrs, and similar to the one given by 
IHardee et all (ITOOl for 3C 449. 

The helical jet structure and a ppearanc e is assumed 
to b e due to KH instab ility as in iHarded (|2000l 120031 ) 
and IHardee et al.1 (|2005[ ). In these works, it is assumed 
that the velocity shear layer is very thin compared to the 
jet radius and wavelength (vortex sheet approximation). 
We note that long wavelengths should be insensitive to 
a modest shear layer so a vortex sheet long wavelength 
analytical approximation would be good enough to de- 
scribe the behavior at long wavelen gths. An alternative 
interpretation to the one in Peruc ho fc Lobanovl (|2007l 
1201 If ) is that the longest wavelength of X'^ ~ 80 mas 
is in the long wavelength and low frequency limit of 
the KH instability spectrum: u <<< uj*, with w* be- 
ing the frequency of the maximum growth rate. Also, 
we might assume that the intermediate wavelength that 



appears to increase from about X'^ ^ 10—12 mas 
{xob ~ — 5 mas) to X°^ ^ 20 — 24 mas {xob 5 — 15 mas) 
to A^'' -- 40 -48 mas {xob ^ 15-35 mas) at the differ- 
ent apparent distances along the jet are each operating 
at about the local resonant wavelength and frequency, 
i.e., a; ~ w* oc aex/rj. Finally, the shortest wavelengths 
X'^ ^ 4.5 mas (xob ~ 5 — 15 mas) and X'^ ~ 7 — 8 mas 
{xob ~ 20 — 30 mas) might be assumed to be a body or 
shear layer mode operating at the appropriate resonant 
frequency and wavelength for that mode. This would 
require the jet having diff erent parameters from those 
obtained in iLobanov et al.l (|l998). 

As we see, different interpretations remain possible: 
vortex-sheet approximation for the longest wavelengths 
or changes in the jet width or in the viewing angle. In 
addition, there are still uncertainties such as the identifi- 
cation in lPerucho fc Lobanovl (|20111 ) of the longest wave- 
length with a first body helical mode, instead of the more 
disruptive, surface mode at the small frequency limit. 
Overall, any results of this kind are subject to sources 
of error associated with the measure of the jet radius, 
the thickness of the shear-layer, the viewing angle or the 
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parameters used for the jet. 

6. SUPERLUMINAL TRANSVERSAL MOTION AND 
NATURE OF THE CORE 

There exist different examples of jets showing apparent 
superluminal tran sversal motion (see, e.g., Agudo et al. 
[20071 [20T0I [20T1I the cases of NRAO 150 and OJ 287). 
We also observe this kind of motion for the ridge-line 
of 0836+710 at 43 GHz. Within the context of the 
present work, there is a possible explanation for this 
phenomenon. We take into account the fact that the 
observations at high frequencies could be showing only 
the high pressure core of the jet along the ridge-line that 
we measure at lower frequencies (see Sections [2] and 14. 1|) . 
This means that we may not observe the whole of the 
jet, but only the region around the high pressure maxi- 
mum. As a result the radio core at 43 GHz would itself 
oscillate transversely and axially over time and not be at 
a fixed position. If the radio core is associated with a re- 
confinement shock ((Marscher et al.|[2008 ). the presumed 
oscillations could well be the triggering point of the ob- 
served growing instabilities, except for the longest one, 
which could be driven by precession at the jet formation 
region. 

Superluminal transversal velocities are possible for 
purely axial motion, a relatively large value for Ay in 
time interval Ai, under some conditions (Perucho et al., 
in preparation). This very special situation may only 
occur for very small viewing angles when Vw is close to 
the speed of light and we deal with high-frequency, short 
wavelength perturbations. There is a second possibility 
that we explore in the next paragraphs. 

6.1. Transversal shift of the apparent core 

One other possible reason for observed superluminal 
transversal motion lies in the mistaken assumption that 
the radio core remains at a fixed (xob, Hob) position. Since 
in our interpretation the radio core is taken as coincident 
with a pressure maximum inside the jet, we should ex- 
pect an oscillation to take place in the radio core loca- 
tion, as the pressure maximum oscillates transversely in 
the plane of the sky and also oscillates axially with rate 
and amount depending on the properties of the helical 
wave (see Eq. 2]) . The oscillation should be small (much 
less than the jet radius) if only a result of an instabil- 
ity amplitude growth but other triggering mechanisms 
can be imagined that could lead to immediate helical 
amplitudes on the order of the jet radius. Therefore, 
when aligning the ridge-line (or any reference point) us- 
ing the jet core as a fixed point at two different epochs, 
it is possible that the whole ridge-line or positions of fit- 
ted components are displaced by the same small amount 
as the core has moved in the transversal and axial di- 
rection. The effect of such artificial displacement will 
clearly affect any measure of velocities in the transversal 
direction. This oscillation has already been reported, 
using phase-referencing , in at least one object (M81, 
IMartf-Vidal"eral][20Tll) . 

For the jet in 0836+710, 1 mas ~ 8.4 pc (see Section[T]), 
and a measured angular motion of only 0.012 mas/yr is 
equivalent to the speed of light. Thus, any small dis- 
placement of the ridge-line anchored at the core that is 
not corrected for possible core motion leads to incorrect 
transversal displacement at any point along the jet and 



can result in bogus measures of superluminal velocities 
in the transverse direction. 

We have generated an artificial mode pattern of a he- 
lical wave with an intrinsic wavelength of A = 2 i?^ , an 
initial amplitude of — 0.2 Rj, and a wave speed of 
0.96 c in order to investigate the role of core shift in the 
measurement of transversal velocities. Figure [T7] shows 
the resulting ridge-lines at two epochs separated by two 
years (left panel) , and the resulting transversal velocities 
in units of the speed of light (right panel). Transversal 
velocities have been obtained by subtracting the vertical 
coordinate of the helix at the different times and dividing 
the result by the time interval, as we did for the obser- 
vational data (see Fig. [15]). With these parameters, the 
core moves transversally at an observed 0.006 mas/yr. 
In this case, the intrinsic displacement of the core is re- 
spected, and we observe that all the velocities are sub- 
luminal with equal positive and negative maximum val- 
ues. The displacement of the ridge-line along the axis is 
a consequence of the projection of the axial coordinate, 
Zh, and that perpendicular to the plane of the sky, Xh, 
(see Eq.[l]). 

In Figure [TSj we re-plot the ridge-lines shown in Fig- 
ure[T7|but now with the cores at the two epochs forced to 
the same position at the origin. This artificial displace- 
ment of the core results in apparent superluminal motion. 
We have investigated the dependence of this effect with 
viewing angle, and axial wave velocity and amplitude. 
The effect does not depend on the viewing angle, as we 
measure transversal motions, but the transversal velocity 
does grow with wave amplitude and wave velocity. Note 
that positive and negative maximum transversal veloci- 
ties are not the same and Figure [18] shows superluminal 
velocities very similar to those found between the 2008 
and 2009 epochs shown in Figure [T^] The similarities be- 
tween the transversal velocities shown in Figure [TS] and 
between the 2008 and 2009 epochs in Figure [TS] provide 
indirect evidence for radio core displacement accompa- 
nying a moving helical wave pattern in 0836+710. A 
direct observational test for core motion would consist of 
phase-referencing observations of near by sources contain- 
ing o bserved helical patterns (see, e.g.. lMartf-Vidal et al.l 
I2OIII) . In these cases small lateral angular displacements 
of the core, up to 0.012 mas/yr if the wave-velocity is 
close to the speed of light, might be directly detectable. 

Core motion manages to explain the transversal veloc- 
ities found in this work in an easy way and within the 
framework of our modeling of helical patterns. We note 
that the observed effect only requires very small displace- 
ments of the ridge-line in the core, well within the errors 
of our measurements. The caveats are that the wave ve- 
locity needs to be close to the speed of light and that the 
initial amplitude of the wave has to be a relatively large 
fraction of the wavelength, 10% in our example. How- 
ever, from stability analysis we know that, at least in the 
case of KH modes, short wavelengths tend to move fast 
and can produce high pressures near to the jet surface. 
Thus, it should not be surprising to find such high pat- 
tern velocities when observing the spine of jets at high 
frequencies and close to the central engine. The proposed 
observational test can be used to verify this hypothesis. 

7. SUMMARY 
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Fig. 17. — Absolute ridge-line position and motion at 3° viewing angle without core position alignment. The left panel shows the ridge- 
lines from the first epoch (solid line) and second epoch (dotted line) separated by a time interval of 2.11 years. The right panel shows the 
resulting transversal velocity versus the observed (projected) axial position. 





Fig. 18. — Ridge-line motion with at 3° viewing angle alignment, the ridge-line from the first epoch is the solid line and second epoch is 
the dotted line. The time interval between them is 2.11 years. On the right, the transversal velocity versus the observed (projected) axial 
position, obtained from the transversal displacement measured at a fixed distance to the core. 



In this work we have presented evidence for the ob- 
served hehcal structures in the jet of 0836-1-710 being 
real, and not artifacts produced by the observational sys- 
tematic errors or to insufficient uv-coverage. The evi- 
dence comes from the same wavelengths being found at 
all epochs, frequencies and observing arrays. We have 
shown that shorter wavelengths develop on top of the 
longer ones and, if their positions along the jet are con- 
sidered to follow a wave structure, the resulting location 
has to be the result of the superposition of the different 
waves and will also depend on the dominant wave at a 
given position and time. 

The fact that shorter wavelengths develop on top of 
longer wavelengths provides a natural way to explain the 
misalignment between parsec and kiloparsec scale heli- 
cal jets seen at small angles to the line of sight and we 
conclude that at the highest frequencies we may be see- 
ing a small region of the jet flow concentrated around 
the pressure maxima that we associated with the inten- 
sity maximum ridge-line. This identification is bolstered 
by 15 GHz VLBA observations, in which we are able to 
resolve the jet, that show the position of the ridge-line. 
This identification is bolstered by 15 GHz VLBA obser- 



vations, in which we are able to resolve the jet, that show 
the position of the ridge-line as not coincident with the 
center of jet. The opening angle of the jet is found to be 
the same at all frequencies, with no apparent frequency 
stratification. 

In the hght of our results, we reinterpret the results 
of IPerucho fc Lobanovl (j2007f ). correcting a mistake in 
the calculations and in the possible identification of the 
triggered instability modes. We confirm the general con- 
clusions obtained in that work to the extent that our 
present study allows but we show that alternative ex- 
planations are still possible. Future work will include a 
detailed discussion on this topic. 

We have pointed out that measuring the velocities of 
the waves is crucial to connecting observed wavelengths 
to the intrinsic wavelengths responsible for the observed 
structures. We have tried to measure the wave veloci- 
ties but the errors associated with the resolution at the 
different frequencies makes this impractical. 

We have also presented a method for modeling helical 
patterns propagating at relativistic speeds in extragalac- 
tic jets. Future work should include the application of 
this method to the observed structures. Finally, we have 
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proposed a possible explanation for moderate transver- 
sal superluminal motions with important implications re- 
garding the nature of the jet-core. In particular, we have 
shown that the radio core at any frequency could be un- 
dergoing transversal motion, implying its possible asso- 
ciation with wave motion. Moreover, the wave-pattern 
obtained for the transversal velocity in the 0836-1-710 jet 
indicates that moving waves are present in the jet. 

Large angular resolution achieved by future space- 
VLBI projects, combined with high precision VLBI as- 
trometry, may help to improve the resolution and make 
it possible to accurately determine the velocities of the 
structures in helical jets in nearby AGN. 
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